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ABSTRACT 


This  study  includes  an  investigation  of  analytical  models  for  the 
description  of  flow  and  bed  characteristics  in  a  bend  of  the  Sacramento  River, 
California;  and  a  study  of  a  proposed  new  technique  for  bank  protection. 

Existing  models  for  the  prediction  of  the  steady-state  transverse  bed 
profile  in  channel  bends  are  reviewed;  and  alternative  models  are  developed 
and  validated.  The  new  models  predict  the  transverse  bed  slope  to  vary 
linearly  with  d/r  (the  ratio  of  depth,  d,  to  radius  of  curvature,  r)  and 
almost  linearly  with  the  sediment-particle  Froude  number.  Equations  for  the 
transverse  variation  of  depth-averaged  velocity  and  mean-grain  size  also  are 
developed,  which,  together  with  the  equations  for  the  transverse  variation  of 
depth,  can  provide  a  guide  to  the  depth  of  bed  erosion  and  depth-averaged 
velocity  near  the  outer  banks  of  river  bends. 

It  is  shown,  theoretically  and  by  a  physical  model,  that  short,  vertical, 
submerged  vanes  installed  at  incidence  to  the  channel  axis  in  the  outer  half 
of  a  river-bend  channel  significantly  reduce  the  secondary  currents  and  the 
attendant  undermining  and  high-velocity  attack  of  the  outer  bank.  The  effect 
of  the  vanes  on  the  secondary  flow  is  estimated  by  a  simple  torque  calculation 
using  the  Kutta-Joukowski  theorem.  A  design  relation  for  the  vane  spacing  is 
derived  by  equating  the  torque,  about  the  channel  centroid,  produced  by  the 
flow  curvature  to  that  resulting  from  the  lateral  force  exerted  on  the 
vanes.  The  relation  is  verified  in  an  idealized,  physical  model  of  the 
mentioned  bend  of  the  Sacramento  River,  California. 
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ANALYSIS  OF  SACRAMENTO  RIVER  BEND  FLOWS,  AND 
DEVELOPMENT  OF  A  NEW  METHOD  FOR  BANK  PROTECTION 


I.  INTRODUCTION 


This  report  presents  the  results  of:  (1)  an  investigation  of  analytical 
models  for  the  description  of  flow  and  bed  characteristics  in  a  bend  of  the 
Sacramento  River,  California;  and  (2)  a  study  of  a  proposed  new  technique  for 
bank  protection.  The  investigation  was  commissioned  by  the  U.S.  Army  Corps  of 
Engineers,  Sacramento  District,  CA  (Contract  No.  DACW05-80-C-0083) ,  and  has 
been  carried  out  in  accordance  with  the  proposal  submitted  to  the  U.S.  Army 
Corps  of  Engineers  by  IIHR  in  December  1979.  The  construction  of  and  testing 
in  a  physical  model  of  mentioned  bend  of  the  Sacramento  River  was  commissioned 
by  the  U.S.  Army  Corps  of  Engineers,  Waterways  Experiment  Station,  MI 
(Contract  No.  DACW39-80-C-0129). 

The  bend  selected  for  this  study  is  located  between  Sacramento  River 
miles  187.5  and  189.5,  which  is  just  above  Ord  Ferry.  The  data  were  obtained 
by  the  U.S.  Geological  Survey  in  April  and  May  1979,  and  March  1980. 

Bank  erosion  is  believed  to  be  caused  primarily  by  flow-induced  boundary- 
shear  stresses.  Because  the  depth-averaged  velocity  tends  to  increase  toward 
the  outer  bank  of  a  river-bend  cross  section,  these  shear  stresses  are  higher 
at  the  concave  bank  than  in  the  central  and  convex-bank  portions.  The  higher 
velocities  at  the  outer  bank,  in  turn,  are  due  to  greater  water  depth  there. 
In  an  attempt  to  provide  a  means  of  predicting  potential  areas  of  bank 
erosion,  this  study  has  investigated  several  mathematical  models  (and 
developed  new  ones)  for  predicting  the  near-bank  velocity.  Since  the  velocity 
distribution  is  intimately  related  to  the  depth  distribution,  much  emphasis 
has  been  put  on  finding  reliable  bed-topography  and  depth  predictors.  The 
effectiveness  of  short,  vertical,  submerged  vanes  as  a  means  of  bank 
protection  also  has  been  evaluated,  both  theoretically  and  in  a  physical  model 
of  the  river  bend. 
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II.  VALIDATION  OF  MATHEMATICAL  MODEL 


A.  Summary  of  Existing  Theories 

Among  the  available  theories  describing  the  flow-bed  interaction  in 
alluvial-channel  bends  are  those  by  van  Bendegom  (1947),  NEDECO  (1959),  Yen 
(1967),  El-Khudairy  (1970),  Ikeda  (1974),  Allen  (1970),  Engelund  (1974a), 
Kikkawa,  Ikeda,  and  Kitagawa  (1976),  Zimmermann  and  Kennedy  (1978),  and  Falcon 
(1979),  All  of  these  theories  are  for  bends  with  fixed  banks  and  constant 
discharges.  They  are  based  on  the  assumption  that  the  mobile  bed  will 
eventually  assume  a  lateral  slope  such  that  the  net  force  on  a  moving  grain  at 
the  bed  will  be  zero  in  the  radial  direction.  The  equilibrium  lateral  bed 
slope  calculated  is  that  for  which  the  radial-plane  components,  parallel  to 
the  bed,  of  the  fluid  drag,  Fp,  and  the  particle  submerged  weight,  Ws,  are 
equal  (Figure  1): 


Fp  =  Ws  sin  0  (1) 

in  which  sin'0  =  radial  (transverse  slope  of  the  bed  surface.  It  is 
principally  in  the  calculation  of  Fp  that  the  theories  differ. 

van  Bendegom  (1947)  determined  Fr  by  balancing  the  forces  acting  on  an 
elemental  fluid  volume  traveling  at  steady  state  through  the  bend.  He 
considered  the  forces  due  to  the  longitudinal  and  transverse  slopes  of  the 
water  surface,  and  the  centrifugal  force.  By  assuming  that:  (1)  the  shear 
stress  behaves  according  to  Prandtl’s  mixing  length  hypothesis;  (2)  the 
direction  of  the  shear  stress  corresponds  to  that  of  the  velocity  vector;  (3) 
the  radial  and  vertical  components  of  the  velocity  vector  are  small  compared 
with  the  longitudinal  component;  and  (4)  the  longitudinal  velocity  profile  is 
described  by  the  Manning-Strickler  power  law  he  obtained 

sin  6  =  K  j  (2) 

in  which  d  =  local  depth  of  flow;  r  =  local  radius  of  channel  curvature;  and  K 
0.0585  F^  ,jS  panicle  Froude  number  defined  as 


F 


D 


V 

S  gD  (Ps-P  )/p 


(3) 
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in  which  p  and  pg  =  density  of  the  fluid  and  sediment,  respectively;  7  = 
depth-averaged  mean  velocity;  D  =  mean  particle  diameter;  and  g  =  acceleration 
of  gravity.  van  Bendegom's  results  were  later  used  by  NEDECO  (1959)  with 
sin  3  replaced  by  the  derivative  d(d)/dr  to  yield  a  differential  equation 
which  was  solved  to  obtain  an  expression  for  the  radial  profile  of  the  bed. 

Yen  (1967),  in  his  analysis  ,  related  Fr  in  Eq.  1  to  the  radial  component 
of  the  flow  velocity,  and  obtained  a  relationship  between  sin  3  and  the 
maximum  radial  velocity.  El-Khudairy 's  (1970)  and  Ikeda's  (1974)  analyses  are 
very  similar  to  the  analysis  by  Yen  (1967),  except  that  the  radial  velocities 
used  in  the  calculation  of  Fp  were  computed  from  Rozovskii's  (1961)  velocity- 
distribution  formulas.  Allen  (1970),  who  was  unaware  of  van  Bendegom's 
results,  used  a  force  balance  similar  to  that  of  van  Bendegom  (1947).  He 
obtained  a  relation  between:  (1)  transverse  slope,  sin  3  ;  (2)  deviation 
angle,  6  ,  between  the  direction  of  the  longitudinal  drag,  Fd,  and  the 
direction  of  the  mean  flow;  (3)  bed  shear  stress,  t  ;  and  (4)  mean  particle 
diameter,  D. 

Engelund  (1974a)  and  Kikkawa,  Ikeda,  and  Kitagawa  (1976)  assumed  that  the 
longitudinal  drag  Fd,  on  a  sediment  particle  is  balanced  by  the  friction  force 
on  the  particle  from  the  channel  bed.  While  Engelund  (1974a)  related  Fr  and 
Fd  by  the  deviation  angle,  5,  using  his  (1974b)  experimental  finding  that 
tan  5  =  7d/r,  Kikkawa,  Ikeda,  and  Kitagawa  (1976)  related  Fr  and  Fd  by 
assuming  that  the  ratio  between  them  equals  the  ratio  between  the  radial,  Vr, 
and  longitudinal,  V0,  velocity  components  of  the  flow  at  bed  level.  V0  was 
determined  from  a  relationship  between  Fd  and  V0,  using  a  drag  coefficient;  Vr 
was  determined  by  solving  the  radial  momentum  equation  using  a  simplified 
stream-function  formulation,  constant  eddy  viscosity,  and  a  logarithmic 
profile  for  the  mean  flow  velocity  (modified  by  the  forced-vortex 
distribution).  Both  Engelund  (1974a)  and  Kikkawa  et  al .  (1976)  obtained  Eq. 
2;  however,  K  is  essentially  a  constant  in  Engelund's  analysis,  while  it  is  a 
function  of  Fq  and  the  Darcy-Weisbach  friction  factor,  f,  in  the  formulation 
of  Kikkawa  et  al.  Engelund's  (1974a)  findings  were  used  by  Engelund  (1974a) 
and  Gottlieb  (1976)  to  calculate  the  bed  topography  of  a  curved  channel  with 
given  discharge  and  grain  size.  Engelund  (1976)  later  showed  that  the 
transverse  bed  slope  increases  when  the  suspended  sediment  load  increases. 
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Zimmermann  and  Kennedy  (1978)  related  Fr  in  Eq.  1  directly  to  stream-flow 
characteristics  by  equating  the  moment,  about  the  channel  axis,  of  the  radial 
boundary  shear  stress  (averaged  around  the  wetted  perimeter)  to  the  moment 
produced  by  the  nonuniformity  (over  depth)  of  the  centrifugal  force.  The 
balance  is  a  simplified  expression  for  the  conservation  of  the  flux  of  moment 
of  momentum.  Falcon's  (1979)  approach  is  an  extension  of  Zimmermann  and 
Kennedy's  (1978).  Falcon  assumed  that  Fr  is  determined  from  a  relationship 
equating  the  centrifugal  force  on  a  unit  volume,  pV2/r,  to  the  radial  pressure 
gradient,  3p/3r,  and  the  vertical  gradient  of  the  radial  shear 
stress,  Sx^/Sy: 


p 


r  3  r  3y 


in  which  p  =  pressure;  y  =  vertical  distance  from  river  bed;  and  x  =  radial 

P 

shear  stress.  The  pressure  gradient  was  determined  by  equating  the  integral 
over  depth  of  the  moment  of  3p/3r  with  respect  to  a  longitudinal  axis,  to  the 
integral  over  depth  of  the  moment  of  the  centrifugal  force  about  the  same 
axis.  Using  the  power  law, 

V  (5) 

in  which  n  =  dimensionless  factor  (which  will  be  termed  the  velocity-profile 
exponent),  the  radial  shear  stress,  x^,  is  determined  as  a  function  of  y  by 
integration.  He  then  obtained  the  radial  drag,  Fr,  by  multiplying  the  radial 
shear  stress  at  the  bed  by  the  projected  surface  area  of  a  sediment  particle, 
A: 


F  =  p  A  V2  N  - 
r  r 


(6) 


in  which 


M  ii'* 

~  n(n+2) 


(7) 


According  to  Zimmermann  and  Kennedy  (1978),  n2  equals  the  inverse  of  the 
overall  friction  factor,  f.  In  Falcon's  (1979)  analysis,  n  is  determined  from 
the  longitudinal  slope,  S,  of  the  water  surface,  using  the  Darcy-Weisbach 
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formula  modified  by  Falcon  (1979).  Both  Zimmermann  and  Kennedy  (1978)  and 
Falcon  (1979)  obtained  Eq.  2,  with 


K  =  cNFg  (8) 

In  Zimmermann  and  Kennedy's  analysis,  the  value  of  c  is  one  third  that  in 
Falcon's  analysis,  because  Falcon  neglected  the  moment  contribution  of 
vertical  shear  stresses. 

Table  1  is  a  summary  of  key  relations  in  the  reviewed  models. 

B.  Alternative  Approaches 

1.  Odgaard's  Approach  -  The  preceding  models  are  based  on  the  assumption 
that  the  equilibrium  bed  profile  is  determined  by  particles  in  downstream 
motion  and  radial  static  equilibrium.  Considering  that  many  natural  river 
beds  (including  those  of  the  Sacramento  River)  consist  of  nonuniform 
sediments,  and  often  have  an  armoring  of  coarser  sediments  at  the  bed  surface, 
there  can  be  a  significant  difference  between  the  grain-size  distribution  of 
the  bed-load  discharge  and  that  of  the  bed-surface  material.  The  following 
approach  (Odgaard,  1981)  is  an  attempt  to  relate  the  radial  bed  slope  to  bed- 
surface  characteristics  rather  than  to  bed-load  characteristics.  Since 
armoring  is  believed  to  be  a  process  that  intervenes  to  bring  the  bed  into 
equilibrium  such  that  no  net  sediment  deposition  or  erosion  occurs,  armoring 
is  considered  to  be  an  incipient-motion  problem.  Therefore,  for  an  armored 
bed,  the  slope  may  be  determined  by  sediment  particles  whose  motion  is 
impending  rather  than  by  particles  that  are  in  motion.  This  approach  is 
basically  a  critical-shear-stress  analysis,  somewhat  similar  to  that  of  White 
(1940)  and  Shields  (1936). 

It  is  assumed  that,  in  the  state  of  equilibrium,  the  sediment  particles 
constituting  the  bed  surface  are  statically  stable  in  both  the  radial 
direction  and  the  longitudinal  direction  (the  mean  flow  direction),  but  at 
incipient  motion  in  the  longitudinal  direction.  That  is,  the  size  of  the  bed¬ 
forming  particles  is  such  that  the  particles  are  just  about  to  move  in  the 
longitudinal  direction  by  rolling  about  their  points  of  support.  Figure  2 
shows  the  vertical  longitudinal-plane  force  diagram  for  one  such  particle. 
Equating  moments  of  the  gravity  (Ws),  lift  (F^),  and  drag  (F^)  forces  about 
the  support  point  yields 
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F 


d 


c2  Ws 


(9) 


in  which 


(a^/a2)tan  <j> 1 
c2  “  1  +  Cj  (a^/'agjtan  <|> 1 

a^  and  a2  are  the  distances  shown  in  Figure  2;  c^  =  the  lift/drag 

ratio,  F^/F^  and  <j>'  =  the  angle  of  repose  of  the  submerged  sediment.  Shields 

(1936)  and  White  (1940)  reduced  Eq.  9  to  the  following  equation  for  the 

critical  shear  stress,  t  : 

cr 


t  =  0  (p  -p  )q  D 
cr  ws  cr 


(ID 


in  which  pg  =  density  of  the  sediment  particles;  g  =  acceleration  of 
gravity;  p  =  density  of  water;  Dcr  =  the  diameter  of  the  particle  whose  motion 
is  impending;  and  e  =  a  dimensionless  factor,  which,  according  to  Shields 
(1936),  is  a  function  of  a  critical  boundary  Reynolds  number.  From  Eqs.  9, 
10,  and  11  it  is  seen  that 

(a1/a2)tan  <j> ' 

0-01  1  +  c1(a1/a2)tan  <j>‘ 

in  which  a  -  ratio  of  projected  surface  are  to  volume  for  a  sediment  particle 

divided  by  that  for  a  sphere.  Shields'  experimental  results,  which  appeared 

as  a  plot  of  t  /[(p  -p)g  D  ]  versus  the  critical  boundary  Reynolds  number, 
cr  s  or 

were  approximated  by  a  curve  (the  Shields  curve)  by  Rouse  (1939).  The  curve 
is  included  in  many  books  on  sediment  transport  (e.g.,  Vanoni,  1970). 
Shield's  curve  is  valid  only  for  flat  beds;  therefore,  xcr  in  Eq.  11  is  the 
critical  flat-bed  shear  stress.  This  shear  stress  may  also  be  related  to  a 
critical  mean-flow  velocity,  Ycp,  by  means  of  the  Darcy-Weisbach  friction 
factor,  f'  (flat-bed  friction  factor): 


T 


_  r 

cr  ~  8 


(13) 


Introduction  of  f'  =  8  tc2n'~^  (obtained  by  comparing  the  power  law,  Eq.  5, 
with  the  logarithmic  relation  for  the  velocity  profile;  Zimmermann  and 
Kennedy,  1978)  yields 
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cr 


x  v  2 


(14) 


If  the  grain-shear  stress,  t',  exceeds  Tcp>  sediment  transport  will  occur. 
However,  substantial  evidence  is  available  (Luque,  1974;  Engelund  and  Fredstfe, 
1976)  indicating  that  the  immobile  bed  will  continue  to  be  acted  upon 
by  x  ,  while  the  residual  part,  t'-t  ,  is  carried  as  drag  on  the  moving  bed 
particles  and  only  indirectly  is  transferred  to  the  bed  by  occasional 
encounters.  If  the  mean  velocity  is  V,  for  a  grain  shear  stress  T'  ,  the 
Darcy-Weisbach  relation  reads 


T  - 


.2 


P  V‘ 


(15) 


and  the  relationship  between  xcr  and  V  may  be  written 

2 


cr 


< _  y2 

~  p  V 
n- 


in  which 


n2  _  n-2  il_  _  .2  V2 

n  r;  —  n  — —  —  n 


cr 


cr 


Equating  Eqs.  11  and  16  yields 


y-2  _  9_  _2  p-"p 


T  nf 

ic 


g  o 


cr 


(16) 


(17) 


(18) 


or 


nf  =  K  F0cr  9 


-1/2 


(19) 


in  which  FQcr  is  the  particle  Froude  number  based  on  Dcr;  that  is, 

V 


Dcr 


g  d 


cr 


(20) 


Falcon  (1979)  and  Zimmermann  and  Kennedy  (1978)  in  their  calculation  of  the 
transverse  bed  slope  assumed  that  the  entire  fluid  stress  is  transmitted 
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directly  through  the  bed-load  layer  to  the  immobile  bed.  In  the  present 

alternative  approach,  the  radial  shear-stress  component  acting  on  the 

particles  of  the  bed  surface  is  the  radial  component  of  the  critical  shear 

stress  t  ,  rather  than  t  .  Therefore,  the  radial  force  balance  is 
cr,r  r 


At  =  W  sin  g 
cr,r  s 


(21) 


in  which  Ws  =  (p$-p)  g  V,  V  =  volume  of  sediment  particle.  In  view  of  Eq.  16 

it  is  reasonable  to  assume  that  x  is  the  radial  bed-shear  stress  necessary 

cr,  r 

to  generate  a  torque  that  balances  the  torque  generated  by  a  power-law 
velocity  profile  of  the  form 

n^+l  .  , 

V  =  $)1/nf  V  (22) 

which,  if  7  =  Vcr,  is  the  velocity  profile  above  a  flat  bed  just  before  the 
onset  of  sediment  motion.  According  to  Falcon  (1979),  the  radial  bed-shear 
stress  is  then  given  by 


u  tt2  d 

Tcr,r  “  pNf  V  7 


in  which 


n^+1 

Nf  =  nf (nf+2) 

Substituting  Eq.  23  into  Eq.  21  yields 

sin  B  =  K  ^ 


(23) 


(24) 


(25) 


where 


(26) 


This  relationship  is  similar  in  form  to  Falcon's  relationship  for  K.  However, 
Falcon's  K  is  a  function  of  the  overall  bed  resistance  and  the  size  of  the 
bed-load  particles,  while  the  present  K-value  (Eq.  26)  is  a  function  of  the 
bed  resistance  given  by  Eq.  19  and  the  size  of  the  bed-surface  particles. 
Because  is  nearly  equal  to  n^  and  n^r,  in  turn,  is  proportional  to  FQCr 
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(Eq.  19),  Eq.  26  predicts  a  nearly  linear  variation  of  K  with  Fpcr.  In 
Falcon's  (1979),  Zimmermann  and  Kennedy's  (1978),  and  van  Bendegom's  (1947) 
models,  K  is  proportional  to  Fp  . 

2.  Falcon  and  Kennedy's  Approach  -  Very  recently.  Falcon  and  Kennedy 
(1982)  proposed  another  approach  for  bed-profile  modelling.  They  used  the 
concept  of  bed-layer  thickness  originally  proposed  by  Karim  (1981)  to 
establish  an  alternative  radial  force  balance  to  Eq.  1.  Considering,  instead 
of  the  single  particle  in  Figure  1,  a  volume  of  bed-layer  particles  of 
thickness  y^,  Falcon  and  Kennedy  proposed  the  following  radial  force  balance 

=  (Ps-p)(l-p')y5  9  sin  3  (27) 

in  which  p'  =  bed-layer  porosity.  Substituting  into  this  equation:  (1) 
Karim's  (1981)  formula  for  the  bed-layer  thickness 

*b  ■  °50  V*/V*c  <28> 

where 

V*  =  V/f/8  (29) 


and 

(vp 

v*c  -yp  g  °5o 9  (30) 

and  (2)  Falcon's  (1978)  Eq.  6  using  n  =  1//7  and  rp  =  Fr/A,  Falcon  and 

Kennedy  obtained  Eq.  2  with  K  determined  by 


1  +  /f  P 
1  +  2/f  U 


(31) 


The  term  in  Eq.  31  involving  f  is  nearly  constant;  therefore,  Eq.  31  predicts 
nearly  linear  variation  of  K  with  Fp. 
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C.  Comparison  of  Predictive  Capabilities 

It  is  seen  that  all  of  the  analytical  models  that  relate  the  local 
transverse  bed  slope  to  mean-flow  characteristics  predict  the  local  transverse 
slope  to  be  proportional  to  the  radio  between  depth,  d,  and  radius  of 
curvature,  r: 


d(d)  _  „  d 
dr  K  r 


(32) 


in  which  Sj  =  d(d)/dr  at  the  centerline. 


1.  Comparison  with  Laboratory  Data  -  Initially,  the  performance  of  the 
models  was  tested  by  means  of  Zimmermann  and  Kennedy's  (1978)  experimental 
results.  Figure  3  (a,b,c,d,  and  e)  shows  a  comparison  of  measured  and 
predicted  slopes,  the  ordinate  being  the  measured  values  of  Sj(dc/rc)  and  the 
abscissa  values  determined  by  the  formula  for  K,  using  measured  values  of  d, 
D,  and  V. 

Figure  3(a)  shows  the  performance  of  the  van  Bendegom  model  (1947).  The 
agreement  between  measured  slopes  and  slopes  predicted  by  this  model 
(calibrated  by  a  factor  of  3.55)  is  quite  good  at  smaller  values  of  Sj. 
Systematic  deviation  of  the  measured  values  of  Sj  from  the  predicted  is  seen 
to  occur  at  the  larger  slopes. 

As  seen  in  Figure  3(b),  the  correlation  between  observed  slopes,  and 
slopes  predicted  by  Engelund's  (1974a)  model,  is  fairly  small.  Since  the 
dynamic  friction  coefficient,  y,  is  more  or  less  constant.  Engelund's  model 
almost  always  predicts  the  same  value  of  Sy/(dc/rc). 

To  apply  Kikkawa  et  al.'s  (1976)  model,  the  function  M  (Table  1)  must  be 
evaluated  at  the  centerline  of  the  channel.  Using:  (1)  A  value  of  one  for  the 
forced  vortex  distribution  (which  is  the  value  to  be  expected  at  the 
centerline);  (2)  v*/V  =  /¥7&  =  l/(n/8);  and  (3)  a  sheltering  factor  of  0.592, 
their  K-value  is  reduced  to  Fg  (1. 628-0. 912/n).  The  value  of  n  for  each 
experiment  was  determined  experimentally  by  Zimmermann  and  Kennedy,  using  the 
side-wall  correction  procedure  of  Vanoni  and  Brooks  (Vanoni ,  1976).  As  seen 
in  Figure  3(c),  the  agreement  between  observed  slopes  and  slopes  predicted  by 
Kikkawa,  et  al.'s  model  (calibrated  by  a  factor  of  0.76)  is  very  good.  Figure 
3(d)  compares  the  observed  slopes  and  the  slopes  predicted  by  Zimmermann  and 
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Kennedy's/Falcon's  model.  The  agreement  is  seen  to  be  somewhat  smaller  than 
when  Kikkawa,  et  al.'s  model  is  applied. 

To  apply  Odgaard's  (1981)  model,  n-f  must  be  evaluated  using  Eq.  19.  The 
value  of  Re*  is  found  to  be  approximately  four  in  Zimmermann  and  Kennedy's 
(1978)  Runs  1-1  through  3-12,  and  approximately  nine  for  the  remaining  runs. 
According  to  Shields  curve,  these  values  of  Re*  correspond  to  values  of  9  of 
0.04  and  0.03,  respecti vely.  As  shown  in  Figure  3(e),  the  agreement  between 
observed  slopes  and  slopes  predicted  by  Odgaard's  model  is  good.  Odgaard's 
model  was  applied  without  introduction  of  a  calibration  factor. 

Figure  3(f)  shows  that  Falcon  and  Kennedy's  (1982)  model  also  conforms 
very  well  to  Zimmermann  and  Kennedy's  data.  The  slope  of  the  straight  line 
through  the  points  gives  a  value  of  porosity  p1  of  0.43  if  Shield's  limiting 
value  of  e  =  0.06  is  adopted. 

2.  Conpari son  with  Field  Data  -  The  ultimate  test  of  a  model's 
predictive  capability  is  against  field  data.  The  data  obtained  by  the  U.S. 
Geological  Survey,  Sacramento,  from  the  bend  of  the  Sacramento  River  between 
river  miles  187.5  and  189.5  were  used  to  verify  the  different  models.  Figure 
4  shows  a  schematic  plan  view  of  the  study  reach.  Six  cross  sections  of  the 
bend  were  surveyed  in  April  and  May  1979,  and  again  in  March  1980  to  measure 
channel  geometry  and  flow  characteristics  for  two  different  flow  rates. 
Detailed  velocity  measurements,  bed-material  samples,  and  depth  soundings  were 
taken  across  each  section.  The  direction  and  magnitude  of  the  flow  velocity 
were  measured  over  4  to  10  verticals  across  each  section,  using  Price  and  Ott 
flow  meters  in  the  1979  surveys,  and  a  Marsh-McBi rney  flow  meter  in  the  1980 
surveys.  In  each  vertical  the  velocity,  V,  was  measured  at  1-  to  2-foot  depth 
intervals;  the  recording  time  was  2  to  5  minutes  at  each  point  of  the 
vertical.  Based  on  these  measurements,  a  depth-averaged  velocity,  "V,  was 
determined  by  dividing  the  known  (gaged)  flow  rate  by  the  measured  cross- 
sectional  area.  The  bed  samples  were  taken  by  a  “digging  sampler"  (USBM-54), 
which  samples  to  a  depth  of  2  to  3  inches  below  the  bed  surface.  Figures  A-l 
through  A-10  in  Appendix  show  log-log  plots  of  V/V  versus  5  =  y/d  for  each 
vertical.  The  velocity-profile  exponents,  n,  were  determined  from  these  plots 
by  measuring  the  slope  of  the  best  fitting  line.  Figures  A-ll  through  A-25  in 
Appendix  show  plots  of  local  and  overall  grain-size  distributions.  (Figures 
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A-24  and  A-25  show  composite  grain-size  distributions  from  an  additional  field 
survey  in  July  1980).  The  results  are  summarized  on  the  drawings  of  the  cross 
sections  in  Figures  5  through  14,  and  in  Table  2.  Two  sets  of  n-values  are 
listed  in  Table  2.  The  values  in  Column  10  were  determined  on  the  basis  of 
the  longitudinal  slope,  S,  of  the  water  surface,  using  the  Darcy-Wei sbach 
formula  modified  by  Falcon  (1979).  The  procedure  for  the  determination  of  n 
is  described  by  Falcon  (1979).  The  slope  of  the  water  surface  was  measured 
and  found  to  be  about  0.0003.  Brice  (1977)  also  reports  that  the  average 
slope  of  the  bed  between  river  miles  190  and  185  is  approximately  0.0003.  The 
n-values  in  Column  11  are  cross-sectional  averages  of  values  determined  from 
the  individual  velocity  profiles  as  described  above.  These  n-values  are  seen 
to  be  significantly  higher  than  those  indicated  by  the  drop  of  the  water- 
surface  elevation  through  the  bend.  That  the  n-values  in  Column  11  are  more 
realistic  characteristics  of  the  bed  friction  than  the  values  in  Column  10  is 
also  indicated  in  Figure  15,  which  shows  a  plot  of  radial  velocity  components, 
Vr,  versus  depth  below  the  water  surface.  The  plot  is  based  on  210  velocity 
measurements  in  the  bend  at  high  flow  (discharges  between  24,000  cfs  and 
28,400  cfs).  Measurements  very  close  to  the  banks  were  omitted.  Each  point 
in  the  plot  is  the  average  value  of  the  radial  velocity  components  over  an 
interval  of  10%  of  the  water  depth.  For  example,  the  value  of  6.9  for 
Vrr'c/(Va  da)  at  K  =  0.9  is  the  arithmetic  average  of  28  measurements 
of  Vrr,~/(VadJ  within  the  range  0.85  <  e,  <  0.95.  The  radial  velocity 

components  were  determined  by  multiplying  the  point  velocities  by  sine  of  the 
angle  between  the  depth-averaged  flow  direction  and  the  flow  direction 
measured  at  each  point  in  the  profile.  The  points  in  Figures  15  are  seen  to 
conform  very  well  to  the  curve  shown,  which  is  the  radial  velocity 
distribution  calculated  from  Falcon's  (1979)  theory  using  n  =  8.5.  The 
average  value  of  n  in  Column  11  (omitting  Section  1)  also  is  8.5. 
Consequently,  frictional  characteristics  and  the  secondary  current  in  this 
bend  may  not  be  estimated  reliably  from  the  longitudinal  slope  of  the  water 
surface.  The  results  presented  in  Figure  15  support  Falcon's  (1979)  theory 
and  the  analytical  approach  used  in  this  study. 

Column  9  in  Table  2  shows  the  transverse  bed  slopes,  Sj.  These  were 
measured  at  the  channel  centerline,  but  they  appeared  to  be  practically 
constant  across  most  of  the  sections.  Figure  16  (a,b,c,d,e,  and  f)  shows  a 
comparison  of  measured  and  predicted  slopes.  The  ordinate  is  the  measured 
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value  of  Sj/(da/rc),  and  the  abscissa  is  the  value  calculated  by  the  formula 
for  K,  using  measured  cross-sectional-average  values  of  velocity,  mean 
particle  diameter,  and  channel  slope;  that  is,  Va,  D5Q,  and  S,  respectively. 
As  will  be  seen  in  Figure  16  (c,e,  and  f),  the  models  by  Kikkawa  et  al . 
(calibrated  by  a  factor  of  2.0),  Odgaard  (calibrated  by  a  factor  of  2.43),  and 
Falcon  and  Kennedy  (calibrated  by  a  factor  of  2.7)  give  good  agreement  between 
measured  and  predicted  slopes.  The  agreement  is  also  good  at  smaller  values 
of  Sj  when  Zimmermann  and  Kennedy's  (calibrated  by  a  factor  of  5.1),  and 
Falcon's  model  (calibrated  by  a  factor  of  1.7),  is  used,  as  will  be  seen  in 
Figure  16(d),  and  when  van  Bendegom's  or  Engelund's  model  is  used,  as  will  be 
seen  in  Figures  16(a)  and  16(b),  respectively;  systematic  deviation  of  the 
measured  values  of  Sj  from  the  predicted  is  seen  to  occur  at  the  larger 
slopes.  Although  Falcon  (1979)  recommends  using  D95  for  the  calculation  of 
Fp,  Falcon's  model  was  found  to  give  better  agreement  between  measured  and 
predicted  slopes  when  D50  was  used.  In  Odgaard's  model,  Nfa  (the  cross- 
sectional  average  value  of  Nf)  was  calculated  by  Eqs.  24  and  19, 
using  e  =  0.06,  <  -  0.4,  and  the  measured  values  of  Fp  for  Fpcr. 

So  far,  the  comparison  of  the  models  has  been  made  on  the  basis  of  their 
capability  of  predicting  the  transverse  bed  slope  at  the  centerline  of  the 
cross  sections.  A  strict  evaluation  of  the  models  should  also  include  a 
comparison  of  the  complete  solutions  in  the  form  recommended  by  the  respective 
authors.  Engelund  (1974),  Kikkawa,  Ikeda,  and  Kitagawa  (1976),  and  Falcon 
(1979)  recommend  using  their  integrated  solutions,  while  Zimmermann  and 
Kennedy  (1978)  and  Odgaard  (1981)  recommend  using  constant  slopes  across  the 
sections.  However,  it  is  easily  verified  that  the  trend  seen  in  Figure  16 
remains  unchanged  whether  integrated  solutions  or  constant  slopes  are  used 
with  the  present  data.  In  all  models,  the  key  parameter  is  the  grain  size, 
whether  it  be  that  of  the  bed  load  or  that  of  the  bed  surface.  As  will  be 
seen  in  Table  2,  the  variation  of  D50  throughout  the  bend  is  very 
significant.  No  technique  or  model  is  yet  available  that  can  simulate  a 
variation  like  that  in  a  given  bend  from  the  flow  and  channel  data.  Instead, 
a  rough  estimate  of  D50  may  be  obtained  on  the  basis  of  a  combination  of  a 
Manning-Strickler  type  of  equation  and  Eq.  19.  Figure  17  shows  a  plot  of  n^a 
(cross-sectional  average  value  of  nf)  versus  relative  depth,  da/D5p.  nf  was 
calculated  by  Eq.  19  using  k  =  0.4,  e  =  0.06,  and  measured  values  of  V  and 
05Q.  It  is  seen  that  a  straight  line  of  slope  1:4.6  gives  a  reasonably 
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accurate  description  of  the  correlation  between  n^a  and  da/D5Q.  The  line  is 
given  by  the  equation. 


nfa  ’  !•«  (da/D50)0'22  (34) 

Substituting  Eq.  34  into  Eq.  19  yields 

D50  =  0.76  da  F3-5  (35) 

in  which  F  is  the  Froude  number,  defined  as 

V 

F  =  (36) 

/gda 

In  Figure  18,  the  measured  correlation  between  relative  roughness  and  Froude 
number  is  compared  with  that  given  by  Eq.  35.  The  agreement  is  seen  to  be 
satisfactory.  Note,  however,  that  Eq.  35  is  based  on  only  data  from  the 
Sacramento  River  study  bend. 

The  procedure  for  application  of  Odgaard's  (1981)  model  is  summarized  as 
follows:  Given  the  cross-sectional  average  values  of  depth,  da,  and  velocity, 
Va,  the  cross-sectional  average  value  of  D50  may  be  estimated  using  Eq.  35. 
(Eq.  35  is  tentative  and  may  apply  only  to  the  Sacramento  River  data).  The 
cross-sectional  average  values  of  n^  and  may  then  be  determined  by  Eqs.  19 
and  24,  respectively.  The  transverse  bed  slope,  sin  3,  is  then  obtained  by 
substituting  n^a,  Nfa,  and  rc  into  Eqs.  25  and  26.  Based  on  the  Sacramento 
River  bend  data,  K  in  Eq.  25  should  be  multiplied  by  a  calibration  factor  of 
2.43. 


The  differential  equation  for  the  transverse  bed  slope,  which  is  obtained 
by  replacing  sin  3  in  Eq.  25  by  d(d)/dr,  is  not  readily  solved  since  both  n^ 
and  FDcr  are  functions  of  mean  velocity  and  grain  size,  both  of  which  vary  in 
a  yet  unknown  manner  across  the  channel  sections.  However,  field  data 
indicate  that  a  reasonable  first  approximation  may  be  obtained  by  considering 
sin  3  to  be  constant  across  a  section.  In  that  case,  the  transverse  depth 
variation  is  described  by  the  equation 


(37) 
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D.  Velocity  and  Grain-Size  Predictor 

It  is  reasonable  to  assume  that:  (1)  the  local  depth-averaged  velocity, 

V,  is  related  to  the  local  longitudinal  slope,  S,  by  the  Darcy-Wei sbach 
equation: 


V  =  n/8gSd 


(38) 


and  (2)  S  varies  as 


/ 


S  =  S 

c  r 


(39) 


in  which  Sc,  rc  =  centerline  values  of  S  and  r,  respectively.  Therefore,  the 
variation  of  the  mean  velocity  across  the  channel  can  be  expressed 

77  .1/2  r  1/2 

zr  =  -ir(lr)  <-r>  (40) 

V  c  c 
c 

in  which  "V  ,  d  ,  n  =  centerline  values  of  V,d,  and  n,  respectively.  Thus,  if 
n  is  constant  across  a  given  section,  that  is,  if  n  =  nc,  the  mean  velocity 
will  vary  with  d/r  to  the  power  0.5.  However,  according  to  the  analysis  by 
Anderson  et  al  (1970)  n  cannot  be  considered  to  be  constant.  They  found  that 
the  Manning  number  for  an  alluvial  channel  varies  with  mean  grain  size  to  the 
1/6-power.  Henderson  (1966)  reports  that  a  1/6-power  relationship  also  was 
observed  by  Strickler  in  this  work  on  gravel-bed  streams.  This  implies  that  n 
varies  with  relative  depth,  d/D,  to  the  power  1/6: 


n  'b 


(41) 


and  Eq.  40  is  reduced  to 


v  D  1/6  2/3  r  1/2 

zr  -  (r>  (f )  ( r*) 

Vc  c 


(42) 


As  can  be  seen  in  Figure  19,  Eq.  41  is  supported  by  the  results  from  the 
Sacramento  River  bend  study.  Eq.  42  shows  that  if  D  =  Dc,  that  is,  if  the  bed 
sediment  is  uniform  and  the  grain-size  distribution  is  the  same  across  the 
width  of  the  channel,  the  mean  velocity,  7,  varies  with  depth,  d,  to  the  power 
of  2/3. 
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If  a  large  variety  of  grain  sizes  is  available  in  the  river  bed,  grain 
sorting  seems  to  take  place.  As  is  illustrated  very  clearly  in  Figures  20  and 
21  (from  Bluck,  1971,  and  Bridge  and  Jarvis,  1976),  the  sediment  near  the 
inner  bank  in  the  downstream  parts  of  the  bends  tends  to  be  finer  than  that  in 
the  deeper  parts  along  the  outer  bank.  The  same  trend  was  observed  in  the 
downstream  parts  of  the  Sacramento  River  study  bend.  Figures  5  through  14, 
from  the  Sacramento  River  bend,  show  that  the  mean  grain  size  can  vary  by  a 
factor  of  more  than  ten  across  a  given  section.  The  variation  of  D  across  a 
section  can  be  estimated  analytically  as  follows:  Eqs.  19  and  40  can  be 
combined  to  yield 


H_  /d. 

n  ^d  ' 
c  c 


1/2  r  1/2  D  1/2  8  1/2 

(f)  <r)  <r> 


(43) 


in  which  nfC  and  9c  are  centerline  values  of  n^  and  0,  respecti vely. 
Substituting  Eq.  43  into  Eqs.  25  and  26  with  Nf  =  n/1,  and  using  Eqs.  40  and 
41,  yields 


5/3  r  3/2  D  2/3 

■  <f>  <r>  <r> 

c 


(S-) 

V 


1/2 


(44) 


in  which  Sjc  =  the  centerline  value  of  Sj.  The  data  from  the  Sacramento  River 
bend,  and  from  several  other  bends  with  nonuniform  bed  sediments,  indicate 
that  the  transverse  bed  slope  may  be  considered  to  be  constant  across  the 
major  part  of  the  cross  section.  From  Eq.  44  it  is  seen  that  in  order  for 
this  to  be  the  case,  the  grain-size  variation  is  governed  by 


D  d  5/2  rc  9/4  0  3/4 

jr  -  <f>  &  (§-> 


(45) 


Substituting  Eq.  45  into  Eq.  42  then  yields 


1/4  r  1/8  er  1/8 
(/)  (jr) 


(46) 


c 

The  Shields1  parameter,  9,  cannot  be  assumed  to  be  a  constant  across  a 
section.  The  data  delineating  Shields'  curve  were  obtained  from  experiments 
in  flumes  with  fully  developed  turbulence  and  artificially  flattened  beds  of 
uniform,  noncohesive  sediments.  The  9-values  were  determined  from  the  values 
of  the  shear  stress  for  zero  sediment  discharge  obtained  by  extrapolating  a 
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graph  of  observed  sediment  discharge  versus  shear  stress.  In  a  river,  the 
flow  conditions  are  very  different  from  those  in  the  flume  experiments. 
Ripples  and  dunes,  formed  as  a  result  of  sediment  movement,  cause  the  value 
of  e  to  be  different  from  that  given  by  the  Shields'  curve  (Shields,  1936; 
White,  1940;  Einstein  and  Barbarossa,  1952;  and  Sundborg,  1956),  and  so  does 
nonuniformity  of  the  bed  sediment.  The  latter  may  be  inferred  from  the  data 
by  Rakoczi  (1975),  who,  in  three  series  of  experiments,  determined  the 
critical  shear  stress  for  different  fractions  of  a  nonuniform  sediment 
mixture.  Figure  22  shows  a  plot  of  his  measured  critical  shear  stresses, 
nondimensional ized  by  (ps~p)g  D^q,  versus  D^q.  Each  point  is  based  on  five 
tests  with  a  given  grain-size  distribution.  The  results  seem  to  indicate  that 
for  a  nonuniform  (non-armored)  sediment  bed,  the  Shields'  parameter  decreases 
with  increasing  mean-grain  size.  The  trend  in  Figure  22  may  be  explained  in 
part  by  Eq.  12  which  shows  that  e  decreases  with  increasing  angle  of 
repose,  <|>'.  According  to  Simons  and  Senturk  (1977),  <j> '  for  coarse  sediment 
decreases  with  increasing  grain  size  within  the  grain-size  range  from  0.5  mm 
to  3.6  mm.  The  decrease  reported  by  them  may  be  described  approximately  by 
the  power  function 


tan  4> '  =  const,  x  D~m  (47) 

in  which  m  =  0.3  and  D  <  3.6  mm.  Substituting  Eq.  47  into  Eq.  12  and  assuming 
that  c^  «  0  yields 


9  ~  D"m  (48) 

Schoklitsch  is  reported  (Graf,  1971)  to  have  suggested  that  m  =  2/3  for  0  <  3 
mm.  As  will  be  seen  in  Graf  (1971),  his  suggestion  was  based  on  a  quite 
significant  number  of  data.  The  trend  which  is  similar  to  that  of  Shields' 
curve  for  Re*  <  10,  may  be  explained  by  viscous  effects:  At  smaller  values  of 
Dgg  a  larger  portion  of  the  sediment  mixture  is  protected  by  the  viscous 
sublayer  than  at  larger  values  of  D50.  Assuming  that  e  can  be  described  by 
Eq.  48  with  m  =  2/3,  which  is  the  value  suggested  by  Schoklitsch  and  used  for 
the  curve  approximating  the  points  in  Figure  22,  and  substituting  Eq.  48  into 
Eq.  45  yields 
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D  d  5/3  rc  3/2 

8-  -  (8-)  o^> 


(49) 


Substitution  of  Eq.  49  into  Eq.  42  then  yields 

V  d  7/18  rc  1/4 

t  >  (r£) 


(50) 


In  Figures  23,  24,  and  25  observed  values  of  (D/Dc)  obtained  in  the  Sacramento 
River  bend  study  are  plotted  against  distance  r  from  the  center  of  curvature, 
and  compared  with  values  computed  by  Eq.  49  using  the  measured  depth 
variation.  The  agreement  is  reasonable.  Figure  26  shows  a  comparison  of 
observed  values  of  V/Vc  and  values  computed  by  Eq.  50,  using  measured  values 
of  d,  r,  V  ,  and  rc.  The  values  plotted  in  Figure  26  are  from  the  central 
portion  of  the  channel  sections  where  bank  effects  are  small.  In  all  cases, 
the  maximum  velocity  measured  in  the  cross  sections  were  within  the  central 
portion  of  the  sections.  Although  the  scatter  is  large,  as  is  invariably  the 
case  with  river  data,  the  results  plotted  in  Figure  26  support  Eq.  50. 


Eqs.  49  and  50  are  based  on  data  from  bend  sections  in  which  the 
availability  of  grain  sizes  were  believed  to  be  unlimited.  Armor  layers  had 
formed  in  many  cases.  If  the  range  of  grain  sizes  in  the  bed  sediment  is 
limited,  the  velocity  distribution  may  be  quite  different  from  that  of  Eq. 
50.  Eq.  42  shows  that  if  D  =  Dc,  the  velocity  distribution  is  given  by 


77  I  2  /3  r  1/2 

Y  V  (r  } 

c  L 


(51) 


If  a  conservative  estimate  of  the  near-bank  velocity  is  needed,  this  equation 
should  be  used  rather  than  Eq.  50. 


E.  Numerical  Example  1 

Given  a  river  bend  of  radius  rc  =  2,000  ft.  Width  b  and  longitudinal 
slope  S  of  the  river  are  500  ft  and  0.0003  ft/ft,  respectively  (based  on 
overall  data);  and  its  overall  frictional  characteristics  are  given  by  n  =  3 
(corresponding  to  an  overall  friction  factor  of  f  =  0.1).  Estimate  the 
maximum  depth,  dQ,  and  maximum  depth-averaged  velocity,  Y  ,  at  a  discharge  of 
Q  =  40,000  cfs. 
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Sol ution: 


Cross-sectional  average  values  of  depth,  da,  and  mean  velocity,  Va,  are 
determined  by  the  Darcy-Wei sbach  relation,  Eq.  38,  and  conservation  of  volume, 
Q  =  Vabda: 


d 


a 


E  Q/b 

n/8gS 


2/3 

■] 


=  21.0  ft 


Vfl  =  n  /8gSd  =  3.8  ft/s 


Grain  size  is  the  major  uncertainty.  If  all  grain  sizes  are  available  in 
the  bed  sediment,  armoring  is  likely  to  occur,  and  D5Q  may  be  estimated  by  Eq. 
35  with  F  calculated  from  Eq.  36: 


°50  =  6,4  01111 

Note  that  Eq.  35  is  based  on  only  the  data  from  the  Sacramento  River  study 
bend. 

The  particle  Froude  number  is  calculated  from  the  grain-size,  D50,  using 
Eq.  20  with  Dcr  =  D50;  V  =  Vg;  and  (ps~p)/p  *  1.65: 


and  nf,  Nf,  and  K  are  calculated  from  Eqs.  19,  24,  and  26/31,  respecti vely , 
using  e  =  0.06;  a  =  1.27;  k  =  0.4;  and  p*  =  0.4: 

n^  =  5.88 

Nf  =  0.148 

K  =  3.65/3.31 

Based  on  the  data  from  the  Sacramento  River  study  bend,  K  must  be  multiplied 
by  a  factor  of  2.43/2.7;  and  the  value  of  K  becomes 


K  =  8.9 
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Maximum  depth,  d0,  and  maximum  depth-averaged  velocities,  7  ,  are  then 

estimated  by  Eqs.  37  and  50,  respectively ,  assuming  that  dc  =  da  , 

and  V  =  V,: 
c  a 


d  =  44.4  ft 
o 

Y  =  4.9  ft/s 
o 

The  estimated  bed  profile  is  sketched  in  the  figure  below.  Note  that  the 
effective  width  of  the  channel  is  reduced  slightly  at  this  discharge.  If  the 
reduction  in  width  had  been  more  significant,  the  calculation  procedure  should 
be  repeated  with  the  reduced  width.  In  effect,  the  solution  then  would  be 
obtained  by  iteration. 


DISTANCE  FROM  OUTER  BANK,  IN  FEET 
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III.  NEW  METHOD  FOR  BANK  PROTECTION 


The  recent  development  in  the  analytical  approaches  has  prompted  the  idea 
that  effective  bank  protection  may  be  obtained  by  means  of  deeply  submerged 
vanes  placed  at  strategic  points  in  the  bend.  From  the  results  of  the 
theories  it  appears  reasonable  to  expect  that  the  strength  of  the  secondary 
flow  in  the  bend  can  be  reduced  if  the  vanes  are  placed  at  a  certain  angle  of 
incidence.  The  vanes  will  then  exert  a  torque  on  the  flow  which  will  counter 
the  torque  produced  by  the  flow  (by  the  interaction  of  the  streamline 
curvature  and  horizontal  vorticity  of  the  flow),  and  a  reduction  of  the 
velocities  along  the  outer  bank  of  the  bend  is  expected. 


A.  Analysis  of  Submerged  Vanes 

The  torque  produced  by  the  flow  through  an  angular  increment  d<j>  may  be 
estimated  using  the  analysis  by  Zimmermann  and  Kennedy  (1978).  Calculating 
the  torque,  dTf|ow,  about  the  center  of  the  flow  yields 

r°  dc  V2  dc 

dTfiow =  /  .  /0  p  r“  (y  -  if)  r  df>dydr  (52) 

in  which  dc  =  depth  of  flow  at  channel  centerline;  V  =  point  velocity  at 

distance  y  from  the  bed  and  distance  r  from  the  center  of  curvature;  r0  = 
outermost  channel  radius;  r-j  =  innermost  channel  radius;  and  p  =  fluid 

density.  The  velocity  V  is  a  yet  unknown  function  of  y  and  r.  Zimmermann  and 
Kennedy  (1978),  however,  assumed  that  the  variation  of  V  with  r  could  be 

neglected  in  the  computation  of  dTf|ow.  They  substituted  the  power  law 

velocity  distribution,  Eq.  5,  into  Eq.  52  and  obtained 


dT 


flow 


=  Bl 
2 


n+1 

n(n+2) 


( r  -  r  • )  d^ 
'  o  i '  c 


V2 


d<f> 


(53) 


This  equation  will  be  used  for  the  present  analysis.  With  reference  to  the 
definition  sketch  in  Fig.  27,  Eq.  53  may  be  re-arranged  to  read 


dTflow  “l^v2^dv 


(54) 


in  which  d¥  =  bdcrc  d<{> ,  the  volume  element  corresponding  to  the  area,  dA, 
hatched  in  Fig.  27(a);  rc  =  (r0  =  rn-)/2;  and  b  =  rQ  -  r^.  The  torque  produced 
by  the  flow  over  a  length  ds  =  rcd<(>  of  the  bend  is 


22 


d 

dTflow  =  2  J(£2TY2b  7^"  ds  (55) 

The  torque  produced  by  a  submerged  vane  may  be  estimated  using  the  Kutta- 
Joukowski  theorem  (see,  for  example,  Daugherty  and  Franzini  1977).  According 
to  this  theorem  the  combination  of  rotational  flow  and  uniform  flow  around  a 
vane  or  foil  of  length  L  and  width  (height  in  the  present  application)  H  gives 
rise  to  a  lift  force,  F£ ,  given  by 


Fz  =  p  HVr  (56) 

in  which  V  =  uniform  flow  velocity  and  r  =  the  circulation  around  the  vane. 
When  the  angle  of  incidence  is  a  ,  the  circulation  is  given  by 


r  =  -rr VL  sin  a 


Combining  Eqs.  56  and  57  yields 


pV2HL 


(57) 


(58) 


in  which 


CL  =  2t r  sin  a  (59) 

Eq.  59  applies  only  to  thin  plates.  Theoretically,  a  may  attain  any  possible 
value;  however,  experience  has  shown  that  a  should  be  less  than  about  20°  in 
order  that  Eq.  59  give  reasonable  estimates  of  cL;  otherwise  flow  separation 
alters  the  pressure  distribution  around  the  foil.  For  a  less  than  about  20°, 
wind-tunnel  experiments  have  yielded  values  of  cL  typically  of  the  order  of 
90%  of  the  values  computed  by  Eq.  59. 

The  torque  about  the  center  of  the  flow  exerted  by  a  vane,  Tyane,  is 

d 

Tvane  Fz  ~  2^  (60) 

in  which  H  =  height  of  the  submerged  vane.  The  height  H  may  realistically  be 
chosen  to  be  of  the  order  of  one  third  the  water  depth;  that  is,  H  =  dc/3 
Hence, 
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^vane  9  a  ^  ^ 

in  which  g  =  ratio  of  actual  to  theoretical  value  of  C|_;  and  V  =  the  flow 
velocity  averaged  over  the  height  of  the  vane.  Averaging  the  velocity  given 
by  Eq.  5  over  the  lower  one  third  of  the  depth  yields 

n+1 

V  =  3(j)n  V  (62) 

The  torque  produced  by  Nv  identical  submerged  vanes,  all  placed  at  the  same 
angle  of  incidence,  is  then 

?  .  (2n+2)/n 

Tvanes  =  sin  «  dc  L  V  Nv  (63) 

In  order  for  the  flow-induced  torque  to  be  balanced  by  the  torque 
produced  by  these  submerged  vanes,  Eq.  63  must  equal  Eq.  55;  that  is 

n+1  b  1  (2n+2)/n 

riXn+57  r  ds  "  2,8  s1n  «  L  Nv  <T>  <64> 

C 

Use  of  typical  data  from  the  Sacramento  River  bend  (b  =  400  ft,  rc  -  2,000  ft, 
and  n  =  8),  and  g  =  0.9,  reduces  Eq.  64  to 

0.0225  ds  =  0.477  L  Ny  sin  a  (65) 

Thus,  if,  for  example,  a  =  20°,  L  =  50  ft,  and  Nv  =  2, 

ds  =  725  ft 

That  is,  for  every  725  ft  along  the  bend,  two  50-foot-long  vanes,  placed  at  an 
angle  of  incidence  of  20°  to  the  local  channel  axis,  should  be  able  to 

counter-balance  the  torque  produced  by  the  flow.  Note  that  this  result  is 

independent  of  the  flow  velocity. 

B.  Physical  Model  Tests 

The  foregoing  theory  has  been  tested  in  a  curved,  recirculating 

sedimentation  flume  at  IIHR,  constructed  with  support  from  the  U.S.  Army  Corps 
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of  Engineers,  Waterways  Experiment  Station,  MI.  The  layout  and  cross-section 
of  the  flume  are  shown  in  Figures  28  and  29.  The  flume  models,  with  some 
idealization,  the  Sacramento  River  bend  at  an  undistorted  scale  of 
approximately  1:48.  The  flume  is  constructed  from  epoxy-coated  plywood. 
Accurately  leveled  2  in.  by  4  in.  plates  on  top  of  the  sidewalls  of  the  flume 
support  a  4  ft  by  8  ft  instrument  carriage.  The  plates  serve  as  vertical 
reference  for  all  elevation  measurements.  Water  enters  the  flume  through  a 
manifold  at  the  upstream  end  of  the  flume  and  leaves  the  flume  at  the 
downstream  end  through  an  adjustable  drop  gate. 

An  approximately  9-in.  thick  layer  of  quartz  sand  was  placed  in  the 
flume.  The  mean  diameter  and  geometric  standard  deviation  of  the  sand  were 
0.30  mm  and  1.45,  respectively.  Water  from  an  underground  sump  was  then 
propelled  through  the  flume  at  a  given  flow  rate  until  the  flow  condition 
became  steady.  Once  a  steady-state  condition  was  obtained  water-surface 
elevations,  and  distributions  of  velocity  and  depth  were  measured.  Velocities 
were  measured  with  a  miniature  propel ler-current  meter,  and  depths  with  a 
sonic  depth  sounder.  Both  current  meter  and  sounder  were  computer  controlled 
(Figures  30  and  31)  and  were  interfaced  directly  with  the  Institute's  HP  1000 
Data  Acquisition  and  Control  System  to  permit  online  processing  of  data. 

C.  Results  of  Physical  Model  Tests 

Numerical  values  presented  in  the  following  are  prototype  values.  The 
model -prototype  scale  ratios  used  are  (Froude  scaling): 

Length:  1:48 
Velocity:  1:6.9 
Time:  1:6.9 
Discharge:  1:15,963 

Figure  32  shows  the  steady-state  bed  topography  measured  in  the  model  at 
a  discharge  (prototype)  of  87,000  cfs,  before  turning  vanes  were  installed. 
The  numbers  by  the  depth  contours  are  depths  in  feet.  The  numbers  along  the 
inside  of  the  flume  are  transverse  bed  slopes  and  section  numbers  (model 
distances,  in  ft,  from  the  manifold,  measured  along  the  right  flume  wall). 
Two  of  the  most  conspicuous  features  of  the  bed  topography  are  the  point  bar 
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along  the  inner  bank  between  Sections  72  and  92;  and  the  relatively  deep  scour 
hole  along  the  outer  (concave)  bank  between  Sections  72  and  100.  The  maximum 
depth  measured  was  41  ft  (at  Section  88).  The  same  features  were  observed  in 
the  Sacramento  River  bend  between  Sections  6  and  3.  As  will  be  seen  in  Figure 
33,  which  shows  the  variations  along  the  Sacramento  River  study  bend  of  water 
level  (WL),  cross-sectional  average  of  bed  level  (BL),  and  level  of  deepest 
points  (SL)  measured  in  the  1980  survey,  a  severe  scour  hole  occurred  between 
Sections  6  and  3  with  the  maximum  depth  increasing  from  17.4  ft  at  Section  6 
to  40.0  ft  at  Section  5  and  34.5  ft  at  Section  4.  At  Section  3  the  maximum 
depth  was  16.5  ft.  The  maximum  transverse  bed  slope  measured  in  the  model  is 
seen  in  Figure  32  to  have  been  0.14  (over  the  point  bar  at  Section  84).  The 
same  slope  was  observed  in  the  Sacramento  River  bend  at  Section  5.  Downstream 
from  the  point  bar  a  very  significant  decrease  of  the  transverse  bed  slope 
occurred.  At  mid-bend  (section  112)  the  transverse  bed  slope  was  reduced  to 
0.02.  This  feature  agrees  with  that  observed  in  the  Sacramento  River  bend 
between  Sections  4  and  3.  Downstream  from  Section  136  in  the  model  the 
transverse  bed  slope  increased  again  to  reach  a  value  of  0.08  at  Sections  144 
through  152.  In  the  prototype,  a  similar  increase  was  not  observed.  This  may 
have  been  due  to  the  heavy  armoring  of  the  prototype  river  bed  between 
Sections  3  and  1,  combined  with  the  increase  of  radius  of  curvature.  Figure 
34  shows  Sections  64,  80,  88,  96,  112,  and  144  in  the  model  before  vanes  were 
installed.  The  n-values  were  determined  from  the  slope  of  the  best  fitting 
line  through  log-log  plots  of  V/V  versus  y/d.  The  average  value  of  n  for  the 
entire  model  bend  was  4.3.  Substituting  n  =  4.3  and  rc/b  -  43/8  into  Eq.  64 
yields 


ds  =  7.3  Ny  L  sin  a  (66) 

Vanes  were  then  installed  in  two  arays  in  the  outer  half  of  the  bend- 
cross  section,  as  shown  in  Figure  35.  Each  vane  was  56  ft  long  and  was  placed 
at  an  angle  of  incidence  of  15°.  At  discharge  87,000  cfs  the  top  of  the  vanes 
was  at  a  depth  of  2/3  the  water  depth,  initially.  The  spacing  between  the 
vanes,  ds  (for  Nv  =  2),  was  determined,  by  Eq.  66,  to  be  212  ft.  The  total 
number  of  vanes  used  was  52.  After  a  period  of  1500  prototype  hours  with  a 
discharge  of  87,000  cfs  the  cross  sections  were  measured  to  be  as '  shown  in 
Figure  36. 
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By  comparing  the  cross  sections  before  (Figure  34)  and  after  the  vanes 
were  installed  (Figure  36)  it  is  obvious  that  the  vanes  had  a  very  significant 
restoring  effect  on  the  bed  topography.  The  point  bar  disappeared  and  so  did 
the  scour  hole.  The  transverse  bed  slope  was  reduced  to  less  than  0.03.  At 
sections  80,  88,  and  96,  where  the  scour  hole  occurred  before  the  vanes  were 
installed,  the  vanes  reduced  the  near-bank  depth  and  velocity  by  approximately 
25%.  The  effect  of  the  vanes  is  even  more  impressive  when  it  is  considered 
that  the  ratio  of  near-bank  to  centerline  velocity  undoubtedly  was 
significantly  higher  in  the  model  than  they  would  have  been  in  the 
prototype.  This  is  due  to  the  smoothness  of  the  flume  wall.  Before  vanes 
were  installed,  a  significant  percentage  of  the  wetted  perimeter  was  flume 
wall  (as  will  be  seen  in  Figure  34)  with  a  much  lesser  friction  than  that  of 
the  sand  bed.  Had  the  roughness  of  the  flume  wall  been  comparable  to  that  of 
the  sand  bed  the  initial  transverse  velocity  gradient  would  have  been  smaller, 
and  fewer  vanes  might  have  accomplished  the  same  bed  changes. 

Very  little  scouring  occurred  around  the  vanes.  Immediately  after  the 
installation  of  the  vanes  some  scouring  was  observed  at  the  upstream  end  of 
each  vane.  However,  as  the  outer  part  of  the  sections  filled  up,  as  a  result 
of  the  vanes,  this  scouring  diminished.  After  a  certain  period  of  time  the 
scouring  was  negligible. 

Tests  also  were  conducted  with  a  smaller  number  of  vanes  (ranging  from  9 
to  52).  It  appeared  that  it  made  little  difference  with  respect  to  velocity 
distributions  whether  the  total  number  of  vanes  was  36  or  52,  or  somewhere 
between,  as  long  as  there  were  two  vanes  per  212  ft  length  (approximately)  in 
the  beginning  of  the  bend  and  at  least  one  vane  per  212  ft  length  in  the 
remaining  part  of  the  bend.  However,  the  transverse  bed  slope  increased 
slightly  as  the  number  of  vanes  was  reduced  toward  36.  Figure  37  shows  the 
steady-state  cross  sections  with  two  rows  of  vanes  in  the  first  part  of  the 
bend  up  to  Section  96  (in  the  configuration  shown  in  Figure  35)  and  one  outer 
row  beyond  Section  96,  giving  a  total  of  36  vanes.  It  is  seen  that  the 
reduction  in  near-bank  velocity  is  about  the  same  as  in  the  case  of  52  anes, 
whereas  the  reduction  in  near-bank  depth  is  slightly  less. 

The  angle  of  incidence,  a  ,  also  was  varied.  For  values  of  a  larger  than 
20°,  flow  separation  occurred  around  30%,  or  more,  of  the  vanes  causing  a 
small,  permanent,  scour  hole  at  these  vanes.  As  the  angle  was  reduced,  the 
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number  of  "scouring"  vanes  was  reduced.  The  optimum  angle,  at  which  no 
significant  scour  occurred,  was  found  to  be  between  15°  and  18°. 

The  overall  effect  of  the  vane  system  on  the  river  flow  was  judged  to  be 
minor.  As  will  be  seen  in  Figure  38,  a  minor,  although  hardly  significant, 
change  of  the  longitudinal  slope  of  the  water  surface  occurred  when  52  vanes 
were  installed  in  the  bend.  The  change  was  less  in  the  case  of  36  vanes.  The 
changes  in  water  depth  also  were  judged  to  be  minor,  and  possibly 
insignificant.  As  will  be  seen  in  Table  3  the  measurements  did  not  disclose 
any  trends  as  to  overall  changes  of  neither  depth  nor  velocity  of  bend  flow  as 
a  result  of  the  vane  systems  tested. 

D.  Numerical  Example  2 

Design  a  vane  system  for  the  bend  in  Example  1,  assuming  that  the  maximum 
flow  is  Q  =  40,000  cfs  during  the  time  from  the  installation  of  the  vanes  to 
full  recovery  of  the  bed  topography. 

The  present  study  did  not  include  optimization  of  the  vane-system 
parameters.  Therefore,  the  configuration  sketched  in  Figure  34  with  L  =  50  ft 
and  Nv  =  2  will  be  adopted.  This  configuration  was  found  in  the  model  study 
to  be  reasonable.  The  value  of  n  to  be  used  in  Eq.  64  should  be  the  value 
based  on  local  flow  characteristics  rather  than  overall  characteristics.  The 
analysis  of  the  Sacramento  River  bend  data  showed  that  the  local  value  of  n 
was  of  the  order  of  8  (whereas  the  overall-value  was  of  the  order  of  3). 
Substituting  the  values 

b  =  500  ft 
rc  =  2,000  ft 
a  =  15° 

0  =  0.9 
L  =  50  ft 
n  =  8 

into  Eq.  64  yields 


ds  =  439  ft 
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That  is,  for  every  439  ft  along  the  bend,  two  50-foot-long  vanes,  placed  at  an 
angle  of  incidence  of  15°  to  the  local  channel  axis,  should  be  able  to  restore 
the  bed  profile.  Installing  the  outer  row  of  vanes  at  a  distance  of  100  ft 
from  the  outer  bank  requires  the  height  H  of  these  vanes  (above  the  mudline) 
to  be 


Houter  “  T~  =  12  ft 

The  maximum  depth-averaged  velocity  at  these  outer  vanes  is  estimated  (by  Eq. 
50)  to  be 


V  = 


(3.8)  (ff) 


7/18 


,2, 000, 1/4 

^2,150 ' 


=4.6  fps 


The  maximum  average  velocity  over  the  height  of  the  vanes  is  estimated  by  Eq. 
62: 


i  9/8 

V  =  (3)  (j)  (4.6)  =  4.0  fps 

The  total  "lift"  force  on  each  vane  is  then  calculated,  by  Eq.  58,  to  be 

F*  =  (7)  (2ir )  (sin  15«  )  (1.94)  (4.0)2  (12)  (50)  =  15,140  lbf 

Installing  the  inner  row  of  vanes  at  a  distance  of  200  ft  from  the  outer  bank 
requires  the  height  of  these  vanes  to  be 

"inner  '  31  ’  8-7  « 

Using  the  same  formulas  as  were  used  for  the  outer  vanes,  the  total  "lift" 
force  on  each  of  the  inner  vanes  is  estimated  at 


=  8,900  lbf 
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IV.  CONCLUSIONS  AND  RECOMMENDATIONS 

The  study  has:  (1)  reviewed  existing  models  for  the  prediction  of  the 
steady-state  transverse  bed  profile  in  channel  bends;  (2)  formulated  and 
validated  alternative  models;  (3)  compared  the  models,  analytically  and  by 
laboratory  and  field  data;  (4)  formulated  and  validated  models  for  the 
prediction  of  the  steady-state  transverse  distributions  of  mean  velocity  and 
grain  size;  and  (5)  evaluated,  theoretically  and  experimentally,  the 
effectiveness  of  vertical,  submerged  vanes  as  a  means  of  bank  protection. 

Eq.  50  has  been  found  to  give  a  reasonably  accurate  relationship  between 
near-bank  velocity  (depth-averaged),  centerline  velocity,  near-bank  depth, 
centerline  depth,  and  the  corresponding  radii  of  curvature.  For  a 
conservative  estimate  of  the  near-bank  velocity,  Eq.  51  is  recommended.  The 
near-bank  depth  may  be  estimated  by  Eq.  37  with  K  being  determined  by  Eq.  26 
or  31.  Based  on  the  data  from  the  Sacramento  River  study  bend,  K  must  be 
multiplied  by  a  factor  of  2.43  if  Eq.  26  is  used;  and  by  a  factor  of  2.7  if 
Eq.  31  is  used.  The  key  parameter  in  all  the  model  is  grain  size.  In  the 
field  surveys  in  the  Sacramento  River  bend  the  mean-grain  size  was  found  to 
vary  significantly  throughout  the  bend.  No  technique  or  model  is  yet 
available  that  can  simulate  or  predict  these  variations  in  a  given  bend,  given 
the  flow  data.  The  grain-size  variation  observed  in  the  field  surveys 
conforms  approximately  to  Eq.  35.  This  equation  may  apply  only  to  the 
Sacramento  River  data.  Due  to  the  importance  of  the  grain  size  for  the 
velocity  and  depth  predictions,  it  is  recommended  that  future  research  be 
focussed  on  this  parameter  and,  in  general,  on  the  bed  processes  of  grain 
sorting  and  armoring. 

The  field  data  from  the  Sacramento  River  study  bend  clearly  show  that  the 
near-bank  velocity,  and  therefore  also  the  intensity  of  the  bank  attack  and 
erosion,  vary  widely  around  the  channel.  For  example,  in  the  1980  survey  the 
depth-averaged  velocity  at  the  outer  bank  at  Section  5  was  measured  to  be  more 
than  100%  higher  than  that  at  Section  3,  2800  feet  downstream  from  Section 
5.  Therefore,  it  is  probably  not  efficient  to  use  the  same  bank-protection 
design  around  a  whole  curve.  Instead,  it  should  be  made  deeper  and  heavier  at 
some  points  (notably  opposite  the  point  bar)  and  tapered  down  to  smaller  depth 
and  thickness  along  other  reaches.  The  analytical  models  developed  in  this 
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study  can  provide  a  guide  to  the  depth  of  bed  erosion  and  depth-averaged 
velocity  near  the  outer  banks,  given  the  discharge,  cross-sectional  area, 
surface  width,  radius  of  curvature,  and  grain  size.  The  models  cannot  predict 
the  formation  of  point  bars. 

A  physical  model  study  is  still  the  most  reliable  means  of  obtaining 
information  about  the  behavior  of  flow  and  sediment  transport  in  a  river 
bend.  A  physical  model  has  been  built  which  models,  with  some  idealization, 
the  Sacramento  River  study  bend  at  an  undistorted  scale  of  approximately 
1:48.  The  model  reproduces  the  main  features  (point  bar,  scour  hole, 
transverse  bed  slope)  very  satisfactorily.  No  mathematical  model  is  available 
that  can  reproduce  these  features  to  an  equally  high  degree  of  accuracy. 
Therefore,  it  is  believed  that  the  best  guide  for  the  placement  of  bank 
protection  is  a  physical  model  study  of  the  bend  in  question. 

The  effectiveness  of  vertical,  submerged  vanes  as  a  means  of  bank 
protection  has  been  evaluated  theoretically  and  in  the  physical  model 
mentioned.  The  studies  showed  that  submerged  vanes,  placed  at  a  certain  angle 
of  incidence  in  the  outer  half  of  the  river-bend  cross-section  have  a  very 
significant,  restoring  effect  on  the  bed  topography  with  a  very  significant 
reduction  of  the  near-bank  velocities  as  a  result.  The  results  indicate  that 
such  vanes  offer  an  effective  (and  environmentally  preferable)  alternative  to 
the  present  bank  protection  measures.  The  vanes  do  not  increase  the  local 
channel  roughness  as  do  other  presently  used  means  for  reduction  of  the  near¬ 
bank  velocity  (for  example,  submerged  dikes);  the  vanes  are  invisible,  and 
they  would  be  far  enough  below  the  surface  of  the  water  that  snags  would  be 
prevented  from  hanging  on  to  them;  and  they  would  appear  to  be  more  economical 
than  traditional  means  of  bank  protection,  especially  bank  armoring. 
Navigation  would  not  be  affected  as  the  vanes  would  be  placed  along  the  outer 
bank  of  the  river  and  the  thalweg  would  move  to  the  central  portion  of  the 
cross  section.  If  either  height  or  total  number  (or  both)  becomes  a  practical 
problem,  lower  and  (or)  fewer  vanes  might  be  combined  with  minimal  slope 
protection  on  the  banks,  such  as  dressing  back  the  banks,  rocking  the  toe,  and 
planting  of  vegatation. 

It  is  recommended  that  the  vane  alternative  be  given  an  early  prototype 
test  after  further  investigation  in  a  model  of  the  specified  bend  in  which  the 
vanes  are  to  be  placed.  The  model  tests  would  be  performed  to  optimize  the 
vane  configuration.  Important  variables  to  analyze  would  be  the  number, 
spacing,  and  length  of  the  vanes  and  their  angle  of  incidence. 
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TABLE  1 

Summary  of  Key  Relations  in  Bed-Topography  Models  Reviewed 


Author  or 

Bed-Friction 

Ratio  of  Local 

Authors 

Parameter 

Bed  Slope  to 
Depth/Radi  us 

Ratio,  K 

(1) 

(2) 

JD _ 

van  Bendegom  (1947) 

C  =  50m1/2/s 

0.0585  Fq2 

Engel und  (1974) 

u  =  tan  $ 

7p 

Kikkawa,  Ikeda, 

M  =  x( 2. 75-4.35  V*/V) f  2 

MF 

Kitagawa  (1976) 

Zimmermann,  Kennedy  (1978) 

N  =  (n+l)/[n(n+2)] 

(a/2)  NFd2 

Falcon  (1979) 

N  =  (n+1 )/  [n(n+2) ] 

(3a/2)  NFd2 

C  =  Chezy's  coefficient 

$  =  dynamic  friction  angle 

fv  =  forced  vortex  distribution 
X  =  sheltering  factor  (  =  0.592) 

V*  =  /r0/p 

a  =  ratio  of  projected  surface  area  to  volume  for  a  sediment 
particle  divided  by  that  for  a  sphere  (=  1.27) 
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Table  2  -  Sacramento 

River  Bend  Flow  Data 

Section 

No. 

Flow  rate,  in 
cubic  feet 
per  second 

Cross-sectional 
area,  in  square 
feet 

Cross-sectional 
average  velocity, 
V  ,  in  feet  per 

Q 

second 

(1) 

(2) 

(3)  

(4) 

1-L 

7,800 

1,900 

4.11 

2-L 

8,300 

3,190 

2.60 

2-H 

24,900 

6,510 

3.82 

3-L 

9,900 

4,340 

2.28 

3-H 

28,400 

7,580 

3.75 

4-L 

9,100 

4,290 

2.12 

4-H 

24,900 

6,370 

3.91 

5-L 

9,100 

4,000 

2.28 

5-H 

24,000 

7,600 

3.15 

6-H _  26.800  6.700  4.00 
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Table  2  (Cont.)  -  Sacramento  River  Bend  Flow  Data 


Section 

No. 

( U 

Cross-sectional 
average  depth, 
d_,  in  feet 

a 

(5) 

Mean  particle 
diameter,  D™, 
in  millimeters 

(6) 

Particle  Froude 
number,  Fn,  based 

on  D50 
(2) 

1-L 

4.4 

21.5 

2.12 

2-L 

5.6 

4.9 

2.81 

2-H 

11.1 

10.8 

2.78 

3-L 

8.1 

6.3 

2.17 

3-H 

12.5 

4-L 

12.2 

1.6 

4.00 

4-H 

23.1 

6.1 

3.78 

5-L 

15.2 

0.7 

6.51 

5-H 

19.7 

- 

- 
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Table  2  (Cont.)  -  Sacramento  River  Bend  Flow  Data 


Section 

No. 

(1) 

Radius  of  curv¬ 
ature,  r  ,  in 
feet  c 

(8)  . 

Transverse  bed 
slope,  Sj 

(9) 

Cross-sectional 
average  velocity 
profile  exponents 
n*^  n** 

no)  on 

1-L 

CO 

0.000 

7.0 

5.8 

2-L 

3,920 

4.0 

6.3 

2-H 

3,920 

0.018 

4.1 

6.2 

3-L 

2,640 

0.020 

2.9 

6.4 

3-H 

2,640 

0.018 

3.8 

5.8 

4-L 

1,800 

0.073 

2.2 

9.6 

4-H 

1,800 

0.135 

2.9 

9.4 

5-L 

1,800 

0.150 

2.3 

10.5 

5-H 

1,800 

0.145 

2.6 

12.5 

6-H 

2,000 

0.041 

4.9 

8.9 

^calculated  from  the  longitudinal  slope  of  the  water  surface 
**) cross- sectional  average  of  values  calculated  from  individual 
velocity  profiles 
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Table  3.  Cross-Sectional  Average  Values  of 
Depth  and  Velocity  in  the  Cases  of  0,  36,  and  52  Vanes. 


Section 

Total  Number 

Cross-Sectional  Average  Value  of 

No. 

of  Vanes 

Depth  (ft) 

Velocity  (fps) 

0 

23.2 

9.8 

64 

36 

24.1 

9.4 

52 

23.6 

9.6 

0 

21.3 

10.6 

80 

36 

21.8 

10.4 

52 

23.9 

9.5 

0 

22.9 

9.9 

88 

36 

24.2 

9.4 

52 

21.4 

10.6 

0 

25.1 

9.0 

96 

36 

22.1 

10.3 

52 

22.3 

10.2 

0 

22.4 

10.1 

112 

36 

21.9 

10.3 

52 

23.6 

9.6 

0 

21.3 

10.6 

144 

36 

22.5 

10.1 

52 

23.4 

9.7 

^Determined  as  the  discharge  per  unit  width  divided  by  the  cross-sectional  average  value 
of  depth. 
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Fig.  1  -  Radial-Plane  Components  of  Forces  on  a  Sediment  Particle 


Figure  2  -  Vertical  Longitudinal -Plane  Components  of  Forces  on  a 
Sediment  Particle. 
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Figure  3  -  Comparison  of  Transverse  Bed 
Slopes  Measured  by  Zimmermann  and  Kennedy 
(1978)  with  Slopes  Computed  by  the  Models 
of  (a)  van  Bendegom  (1947);  (b)  Engelund 
(1974);  (c)  Kikkawa,  Ikeda,  and  Kitagawa 
(1976);  (d)  Zimmermann  and  Kennedy  (1978) 
or  Falcon  (1979);  and  (e)  Odgaard  (1981) 
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Figure  5  -  Cross  Section  No.  1  at  Low  Flow,  7,800  Cubic  Feet  Per  Second 
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Figure  6 


Cross  Section  No. 


2  at  Low  Flow,  8,300  Cubic  Feet  Per  Second 
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Figure  7  -  Cross  Section  No.  3  at  Low  Flow,  9,900  Cubic  Feet  Per  Second 
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Figure  8 


Cross  Section  No.  4  at  Low  Flow,  9,100 
Cubic  Feet  Per  Second 
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Figure  9 


Cross  Section  No.  5  at  Low  Flow,  9,100 
Cubic  Feet  Per  Second 
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Figure  10  -  Cross  Section  No.  2  at  High  Flow, 
24,900  Cubic  Feet  Per  Second 
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Figure  11  -  Cross  Section  No.  3  at  High  Flow, 
28,400  Cubic  Feet  Per  Second 
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Figure  12  -  Cross  Section  No.  4  at  High  Flow, 
24,900  Cubic  Feet  Per  Second 
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Figure  13  -  Cross  Section  No.  5  at  High  Flow,  24,000 
Cubic  Feet  Per  Second 


Figure  14  -  Cross  Section  No.  6  at  High  Flow,  26,800 
Cubic  Feet  Per  Second 
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Figure  15  -  Measured  and  Computed  Radial  Velocity  Distributi 
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Figure  16  -  Comparison  of  Transverse 
Bed  Slopes  Measured  in  Sacramento  River 
Bend  with  Slopes  Computed  by  the  Models 
of  (a)  van  Bendegom  (1947);  (b)  Engelund 
(1974);  (c)  Kikkawa,  Ikeda,  and  Kitagawa 
(1976);  (d)  Zimmermann  and  Kennedy  (1978) 
or  Falcon  (1979);  and  (e)  Odgaard  (1981) 
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Figure  16(f)  -  Comparison  of  Transverse  Bed  Slopes  Measured 
in  the  Sacramento  River  Study  Bend  with  Slopes 
Computed  by  Falcon  and  Kennedy's  (1982)  Model 
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Figure  19  -  Velocity  Profile  Exponent  n  as  Function  of  Relative  Depth,  Sacramento  River 


Row  Direction 


Figure  20  -Distribution  of  Sand  and  Gravel  on  the  Bed  of  River 
Endrick,  Scotland  (Bluck  1971  ). 


Figure  21  -  Distribution  of  Sand  and  Gravel  on  the  Bed  of  River 
Esk,  Scotland  (Bridge  and  Jarvis  1976  ). 
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Fig.  22  -  Shields'  Parameter  as  a  Function  of  Mean-Grain 
Size,  Based  on  Rakoczi's  (1975)  Measurements. 
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Figure  23  -  Observed  and  Predicted  Transverse  Variation  of  Grain 
Size  at  Section  3  of  the  Sacramento  River  Study  Bend. 
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Figure  24  -  Observed  and  Predicted  Transverse  Variation  of  Grain  Size 
at  Section  4  of  the  Sacramento  River  Study  Bend 
(1  ft  =  0.305  m). 
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Figure  25  -  Observed  and  Predicted  Transverse  Variation  of  Grain 
Size  at  Section  5  of  the  Sacramento  River  Study  Bend 
(1  ft  =  0.305  m). 
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039  0  25 

COMPUTED  DIMENSIONLESS  MEAN  VELOCITY,  (d/dc)  (rc/r) 


Figure  26  -  Comparison  of  Measured  and  Computed  Transverse  Variation 

of  Mean  Velocity  across  Bends  of  Desna  River  (Rozovskii  (1961)), 
Kaskaskia  River,  Illinois  (Bhowmik  (1979)),  and  Sacramento 
River,  California. 
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FORCE  PROFILE 


SIDE  VIEW 
SECTION  A-A 


SL 


tvane 

■>  - 

fj  VANE 

SIDE  VIEW 
SECTION  B-B 

(b) 


(a) 

Figure  27  -  Definition  Sketch  for  the  Calculation  of  the 
Torque  Produced  by  (a)  the  Flow  and  (b)  a 
Submerged  Vane 
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Figure  28  -  Layout  of  the  Curved,  Recirculating  Sedimentation 
Flume 


Figure  29  -  Cross  Section  of  the  Curved,  Recirculating  Sedimentation 
Flume 


Figure  30  -  Computer  Operated  Depth  and  Velocity  Measurements 
in  the  Curved  Sedimentation  Flume 


Figure  31  -  Close-up  of  Current  Meter  and  Sonic  Depth  Sounder 
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Figure  33  -  Variations  along  the  Sacramento  River  Bend  of 

Water  Level  (WL),  Cross-Sectional  Average  of  Bed 
Level  (BL),  and  Level  of  Deepest  Point  (SL), 
Measured  in  the  1980  Survey 
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Figure  36 


Steady-State  Cross  Sections  with  52  Vanes 
Installed;  Discharge  =  87,000  cfs 
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Figure  37  -  Steady-State  Cross  Sections  with  36  Vanes 
Installed;  Discharge  =  87,000  cfs 
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Figure  A-2  (Cont)  -  Velocity  Distributions  in  Cross  Section 

No.  2  at  Low  Flow,  8,300  Cubic  Feet  Per 
Second 
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Figure  A-3  -  Velocity  Distributions  in  Cross  Section  No.  3 
at  Low  Flow,  9,900  Cubic  Feet  Per  Second 
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Figure  A-4  -  Velocity  Distributions  in  Cross  Section  No.  4 
at  Low  Flow,  9,100  Cubic  Feet  Per  Second 
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Figure  A-6 


Velocity  Distributions  in  Cross  Section  No.  2 
at  High  Flow,  24,900  Cubic  Feet  Per  Second 
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n  *  logH  (1.18/0.8)  *  5.9 


Figure  A-7  -  Velocity  Distributions  in  Cross  Section  No.  3 
at  High  Flow,  28,400  Cubic  Feet  Per  Second 
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Figure  A-10  -  Vertical  Velocity  Distributions  in  Cross 
Section  No.  6  at  High  Flow,  26,800  Cubic 
Feet  Per  Second 
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Figure  A-11  -  Grain-Size  Distribu¬ 
tions  in  Cross  Section  No.  1  at 
Low  Flow,  7,800  Cubic  Feet  Per 
Second 


Figure  A-l 2  -  Grain-Size  Distributions  in 
Cross  Section  No.  3  at  Low  Flow,  9,900  Cubic 
Feet  Per  Second 
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Figure  A-l 3  -  Grain-Size  Distributions 
in  Cross  Section  No.  2  at  Low  Flow, 
8,300  Cubic  Feet  Per  Second 


Figure  A-14  -  Grain-Size  Distributions  in 
Cross  Section  No.  4  at  Low  Flow,  9,100 
Cubic  Feet  Per  Second 
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Figure  A-15  -  Grain-Size  Distribu-  Figure  A-16  -  Composite  Grain-Size  Distribu¬ 
tions  in  Cross  Section  No.  5  at  tions  at  Low  Flow 
Low  Flow,  9,100  Cubic  Feet  Per 
Second 
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Figure  A-17  -  Composite  Grain-Size  Distribution  for 
all  Sections  at  Low  Flow 
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Figure  A— 1 8  -  Grain-Size  Distributions 
in  Cross  Section  No.  2  at  High  Flow, 
24,900  Cubic  Feet  Per  Second 


Figure  A-19  -  Grain-Size  Distribu¬ 
tions  in  Cross-Section  No.  4  at  High 
Flow,  24,900  Cubic  Feet  Per  Second 


Figure  A-20  -  Grain-Size  Distribu¬ 
tions  in  Cross  Section  No.  5  at  High 
Flow,  24,000  Cubic  Feet  Per  Second 


Figure  A-21  -  Grain-Size  Distributions 
in  Cross-Section  No.  6  at  High  Flow, 
26,000  Cubic  Feet  Per  Second 
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Figure  A-22  -  Composite  Grain-Size 
Distributions  at  High  Flow 


Figure  A-23  -  Composite  Grain-Size 
Distributions  for  all  Sections  at  High 
Flow 


Figure  A-24  -  Composite  Grain -Size 
Distributions,  from  Survey  in 
July,  1980 


Figure  A-25  -  Composite  Grain-Size 
Distribution  for  all  Sections,  from 
Survey  in  July,  1980 


